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Abstract 

We consider the step Riemann problem for the system of equations describing the propagation 
of a coherent light beam in nematic liquid crystals, which is a general system describing nonlinear 
wave propagation in a number of different physical applications. While the equation governing the 
light beam is of defocusing nonlinear Schrodinger equation type, the dispersive shock wave (DSW) 
generated from this initial condition has major differences from the standard DSW solution of the 
defocusing nonlinear Schrodinger equation. In particular, it is found that the DSW has positive 
polarity and generates resonant radiation which propagates ahead of it. Remarkably, the velocity 
of the lead soliton of the DSW is determined by the classical shock velocity. The solution for the 
radiative wavetrain is obtained using the WKB approximation. It is shown that for sufficiently 
small initial jumps the nematic DSW is asymptotically governed by a Korteweg-de Vries equation 
with fifth order dispersion, which explicitly shows the resonance generating the radiation ahead of 
the DSW. The constructed asymptotic theory is shown to be in good agreement with the results 
of direct numerical simulations. 


1 Introduction 

Dispersive shock waves (DSWs), also termed undular bores in fluid mechanics, are generic solutions 
of nonlinear dispersive wave equations, including the Korteweg-de Vries (KdV), nonlinear Schrodinger 
(NLS) and Sine-Gordon equations. A DSW forms due to the dispersive resolution of a discontinuity and 
is the dispersive equivalent of a gas dynamic shock for which a discontinuity is resolved by viscosity [T] . 
A DSW is a non-steady modulated wavetrain which continually expands and has solitary waves at its 
leading edge and linear, small amplitude waves at its trailing edge (for the case of negative dispersion; 
if dispersion is positive then the orientation of the DSW, i.e. the relative position of the linear and 
soliton edges, changes). This modulated wavetrain provides an oscillatory transition between the two 
levels of the initial discontinuity. 

DSWs/undular bores are a common wave form which can be found in a broad array of physical 
systems. The classical undular bore is the tidal bore found in regions of large tidal flows and suitable 
topography, for example the Severn Estuary in England and the Bay of Eundy in Canada. However, 
undular bores arise in a wide range of fluid systems, including the atmosphere, an example being 
morning glory clouds mE], and the semi-diurnal internal tide [4]. They also arise in geophysics 
(magma flow) [71 EJ [5] and Fermi gases |8]. Of particular relevance to the present work, they arise in 
nonlinear optics for a wide range of optical materials, including photorefractive crystals mi El do], 
optical fibres mmsiiiiiiiaiiniiiT], nonlinear thermal optical media [IH11I9], colloidal media I201EI] 
and nematic liquid crystals [2I1122]. 

DSW solutions of nonlinear dispersive wave equations are usually found using Whitham modulation 
theory [231 1211 IB- Whitham modulation theory is a method for analysing slowly varying (modulated) 
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wavetrains and deriving equations for the parameters, mean height, wavenumber, amplitude, etc., of 
such wavetrains. It is equivalent to the method of multiple scales, but much simpler than this to 
implement. When the underlying wavetrain is stable, the modulation equations form a hyperbolic 
system for the wavetrain parameters. It was found that a simple wave solution of the hyperbolic mod¬ 
ulation equations for the KdV equation corresponds to a DSW [25l|26], so that the standard method 
for finding DSW solutions is from the modulation equations for the relevant governing equation. This 
original method due to Gurevich and Pitaevskii [25] and Fornberg and Whitham [26| relies on the 
hyperbolic modulation equations being in Riemann invariant form, which is guaranteed if the govern¬ 
ing equation is integrable with an inverse scattering solution m- However, most equations governing 
DSWs in physical applications are not integrable. This limitation was overcome to a certain extent 
when it was found that the leading, soliton, edge and trailing, small amplitude wave, edge of a DSW 
could be determined without a knowledge of the full Whitham modulation equations |28j . 

In the present work a DSW due to coherent light propagation in a nematic liquid crystal is analysed. 
While the specihc context is light propagation in a nematic liquid crystal, equations similar to those for 
light propagation in this medium also arise for other nonlinear optical media [snilHIllSgllMl [aSllMl [E] , 
in fluid mechanics [35| and in models of quantum gravity [36]. An optical DSW in a nematic liquid 
crystal is found to possess a number of unique features. While the equation governing the optical 
field in a nematic liquid crystal is of defocusing NLS-type 133, the DSW is found to be of positive 
polarity, KdV-type, due to the effect of the nematic optical medium, which has a highly “nonlocal” 
response |391 IdOt 138]. It is further found that the dispersion relation for linear waves is non-convex, 
so that there is a resonance between the DSW and dispersive radiation. This results in a resonant 
wavetrain propagating ahead of the DSW. A similar resonant coupling between a DSW and radiation 
was found for nonlinear optical beam propagation in optical hbres when higher order dispersive terms 
were included in the governing NTS equation to enable such coupling, both without Haim [17] and 
with m loss. The driving mechanism is the resonant coupling with higher order dispersion, which 
can also occur with just a soliton m- The total structure of the Riemann problem solution is then 
found to consist of four distinct regions, (i) an expansion wave linking the initial level behind to an 
intermediate shelf, (ii) a KdV-type DSW on this shelf, (hi) a resonant wavetrain leading the DSW 
and (iv) a front bringing the resonant wavetrain down to the initial level ahead. Asymptotic solutions 
for all these four regions are obtained and compared with full numerical solutions of the governing 
equations, with generally excellent agreement being found. 

The paper is organised as follows. In Section [2] the equations governing light beam propagation in 
a nematic liquid crystal are introduced and related to similar systems of equations in other physical 
contexts. In Section [3| the dispersive-hydrodynamic properties of these nematic equations are analysed 
and it is found that, while the dispersionless limit is described by a hyperbolic system equivalent to 
the shallow water equations, which is consistent with the dispersionless limit of the defocusing NLS 
equation, the linear dispersion relation is non-convex, implying the possibility of the formation of a 
KdV-type DSW in the low frequency region and the generation of high frequency resonant radiation 
by the DSW. This effect is a counterpart of the well known radiating solitary waves in systems with 
higher order dispersion studied previously in many physical contexts, from gravity-capillary waves 
[42] to optical supercontinuum generation (see e.g. [43] and references therein). In Section [4| the hfth 
order KdV equation (also known as the Kawahara equation) is derived from the nematic equations 
under a balance between strong nonlocality and the small amplitude, long wave approximation. The 
coefficient of the fifth order dispersion term is proportional to the nonlocality squared. It is then 
shown numerically that the effect of the nonlocality on the DSW is the generation of a radiative 
wavetrain ahead of the DSW. In contrast to the well studied radiating solitons of the hfth order KdV 
equation, which are intrinsically unsteady, the solitary wave at the leading edge of the radiating DSW 
remains steady due to energy inhux from the rest of the DSW. It can then be well approximated 
by the standard KdV soliton if the higher order dispersive term is sufficiently small. In contrast to 
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previous work [22] it is found that the velocity of the leading edge of the KdV DSW is given by a 
classical shock jump condition, rather than the conservation of Riemann invariants [28|. This suggests 
that the resonant wavetrain acts as an effective viscous loss term for the DSW. In Section 05.31 a 
WKB solution is constructed for the rapidly oscillating, resonant, linear radiative wavetrain in the full 
nematic system under the assumption that the lead solitary wave in the DSW can be approximated 
by a KdV soliton. Section [6] is devoted to comparisons of the constrncted modulation solution with 
fnll numerical solutions of the nematic system. 

2 Nematic equations 

In this paper, we consider the propagation of a polarised, coherent beam of light through the medium 
of a nematic liquid crystal [4f)t 138] . We assume that the electric field of the light is in the x direction 
and that the beam propagates in the z direction. Nematic molecules are elongated molecules, hence 
their name as nematic comes from the Greek word for thread, along which electrons can move freely. 
Hence an electric field, either an external static electric field or the electric field of light, results in 
the nematic molecules becoming dipoles and rotating in the direction of the electric field due to the 
resulting torque in order to minimise the potential energy I1Q1138]. The molecules rotate until the 
elastic forces balance the electrostatic forces. This rotation changes the refractive index of the nematic 
medium. Normally a nematic is a focusing medium, so that the refractive index increases on rotation 
of the molecules. This self-focusing can then balance the diffractive spreading of a light beam, so that 
a bright optical solitary wave, termed a nematicon, can form [441 HOl 138] . However, the addition of 
azo-dyes to the nematic medium changes its structure so that it can become a defocusing medium as 
rotation of the molecules then decreases the refractive index |l5]. In this case, a dark solitary wave, 
a dark nematicon, can form, a dip in a uniform background, rather than the rise from a background 
of a bright nematicon in the focusing case. The added complication of the nematic medium is that 
if the nematic molecules are initially aligned with their axis, termed the director, orthogonal to the 
electric field, the optical Freedericksz threshold exists so that a minimum electric field strength is 
required to overcome the elastic forces of the nematic medium before the molecules can rotate [461 00] . 
To enable nematicons to form at milliwatt power levels so that there is not excessive heating of the 
nematic, which can result in it undergoing a phase transition, an external static electric field is applied 
to pre-tilt the nematic molecules at an angle 9q to the 2 : direction. In the particular case 6q = 7r/4, 
the Freedericksz threshold vanishes [441 fSO] . 

Let us denote the extra rotation from the pre-tilt caused by the electric field of the light beam to 
be 9. Then in the paraxial, slowly varying envelope approximation the system of equations governing 
the propagation of a nonlinear light beam through a defocusing nematic liquid crystal is [3^051001138] 

0 , 

-2\u\^. ( 2 ) 

Here u is the complex valued envelope of the electric field of the light beam. The parameter u, termed 
the nonlocality, measures the elastic response of the nematic and is large, v = 0(100), in experiments 
m- This large value of the nonlocality i' will be found to have a dominant effect on the structure 
of a DSW in a defocusing nematic liquid crystal. The parameter q is proportional to the square of 
the pre-tilting electric field. The electric field equation ([1]) is a nonlinear Schrodinger (NLS)-type 
eqnation, which is coupled to equation ([2]) for the response of the nematic medium. 

The context of the system of equations ([T]) and ([2]) has been explained in detail in terms of the 
nonlinear optics of liquid crystals. However, this system arises in a wide range of applications. In 
nonlinear optics, it arises whenever the response of the optical medium is based on some type of 
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Figure 1: Numerical solution of nematic equations dH) and ([5]) for = 1.0, ui = 0.5 and q = 2. Initial 
condition for |u| (z = 0): red (short dashed) line; initial condition for 6 (z = 0): green (long dash) 
line; numerical solution for \u\ at z = zf. blue (solid) line; numerical solution for 9 at z = zf. pink 
(dotted) line, (a) v = 200, Zf = 600, inset detail of resonant wavetrain, (b) u = 0.1, Zf = 300. 


diffusive phenomenon [29], for example it arises in the optics of nonlinear thermal media 

for example lead glasses [sniEaEsi, and certain photorefractive crystals [34|. A similar system of 

equations arises in simplified models of fluid turbulence |35| and quantum gravity |36j . 

In this paper we consider the Riemann problem for the nematic system ([T|) and ([2|). The electric 
field equation ([1]) will be solved with the initial condition 


u = 


Us, X < 0 
ui, X > 0 


(3) 


at z = 0, with Us > ui so that a DSW is be generated. For consistency, the director equation ([2|) gives 


03 = ^, X<0 

01 = ^, x>0 


(4) 


at z = 0. 

A typical solution of the nematic equations for the step initial condition ([Hj) for large nonlocality 
u is displayed in Figure [iKa). For comparison, the solution for small v is displayed in Figure [iKb), 
noting that the nematic equations ([1|) and Q reduce to the NLS equation in the limit —)• 0. As 
found in previous work [22|, for large values of the nonlocality v the solution does not display the 
typical defocusing NLS DSW structure of Figure [I](b) [37], even though the electric field equation 
([T]) is of defocusing NLS-type. There is a KdV-type DSW in the electric field on the intermediate 
shelf of height U 2 between the initial levels us and ui. Preceding this DSW, there is a relatively 
high frequency wavetrain, with a front which brings it back to the initial level ui. The KdV DSW 
and resonant wavetrain are mirrored in the director response, at a much reduced amplitude, with 
the resonant wavetrain in the director having amplitude The inset in Figure [T]( a) shows the 

details of this resonant wavetrain. In this paper, the complex wave structure seen in Figure [11(a) is 
understood as a radiating DSW. Such radiating DSWs typically arise for nonlinear wave equations 
with higher order dispersion, the model equation being the fifth order KdV equation, or Kawahara 
equation. Although the theory of radiating solitons for the fifth order KdV and similar equations 
is well understood, see e.g. [1H11191I50] and references therein, the counterpart for DSW theory has 
only started to be explored (see the monograph m and references therein and the recent papers 
[la la la isi dz]). In addition to the leading resonant wavetrain, there are also radiative waves on 
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the intermediate shelf on which the DSW sits. These are most likely due to internal resonances within 
the DSW which is a modulated periodic wave with a range of phase and group speeds. 


3 Nematicon dispersive hydrodynamics 


To analyse the Riemann problem ([I])-(IH) it is instructive to introduce the Madelung transformation 


u = V = (l)x 


(5) 


in order to set the nematic equations 


m and ([2]) in the so-called dispersive hydrodynamic form 


dv dv do 
dz^ ^ dx^ dx 
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( 6 ) 

(7) 

( 8 ) 


The above hydrodynamic form highlights the presence of two characteristic spatial scales in the system 
for large v. the long scale and the short scale 0(1), which is consistent with the two distinct 

types of oscillatory structures observed in Fig. [11(a). These distinct structures are characterised by 
differing typical wavelengths and different types of dispersion, which can be understood by analysing 
the linear dispersion relation for this system. 

Linearising the hydrodynamic form of the nematic equations (l6])-(|8|) around the background levels 
p, V and 9 with 

p = p + p, V = V + V, 9 = -+ 9, (9) 

Q 

where \p\ <C p, |u| <C |u| and |0| <C p/q, gives the dispersion relation for right-propagating waves [22] 


uj = kv + 


y/f^ 


y/uW+2q [ 4/3 


vk'^ 2q 


k^ + Ak 


nl/2 


( 10 ) 


We note that since the dispersion relation (ITOl) is obtained not for the original system ([T]) and (|2|), but 
for its dispersive-hydrodynamic representation (I6|)-(l8|), it does not contain the frequency shift 20i 
due to the background carrier wave Y^exp(— 2 i 0 i 2 ;). 

To better understand the dispersive properties of the nematic system given by the dispersion 
relation (|10p we consider its long wave and short wave expansions. Expanding (|10l) in powers of A: <C 1 
and retaining terms up to 0{k^) we have 


UJ ~ k[c + v) — 



4pj 32 V ^ p 16p2 J 


( 11 ) 


where c = ^J2pjq. The expansion dill) requires not just that A: <C 1, but that uk"^ = 0(1) or uk^ <C 1, 
which generally does not hold true even for reasonably small wavenumbers k due to the very large value 
of the nonlocality v. Nevertheless, as we shall see, the expansion ([TT|) captures some key qualitative 
features of the full dispersion relation. Now looking at the short wave asymptotics of (jinj) . we obtain 
that for strong nonlocality v ^ 1 

UJ zz kv+ ^k^ + 0{{k^/v)~^), vk^ ^ 1. (12) 
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Figure 2: Nematic dispersion relation (a) Full dispersion relation cni): red (solid) line, 5th order 
expansion (HU): blue (dashed) line, (b) nematic phase velocity C = oj/k: red (solid) line, 4th order 
expansion velocity: blue (dashed) line. The parameters are p = 1.0, v = 0, u = 200 and q = 2. 


One can see from the expansions m and m that for sufficiently small wavenumbers ujkk < 0, 
while for large wavenumbers LOkk > 0. Thus the full dispersion relation (|10l) is non-convex, which 
has important physical consequences as it implies the possibility of resonance between long and short 
waves and hence the generation of short wave radiation propagating ahead of the DSW. The effect of 
resonant radiation generation by solitary waves in equations with higher order dispersion is well known 
in the context of gravity-capillary waves, see e.g. [m ig Hg [50] and references therein. There is also 
abundant literature on radiating solitary waves in nonlinear optics, see e.g. |52[ l43] and references 
therein. However, the counterpart of this for DSW theory is yet to be developed. A few existing 
notable contributions include numerical investigations of radiating DSWs described in the monograph 
m and the recent papers [m o Esiiini IT] on the effects of higher order dispersion on NTS DSWs 
in the context of nonlinear optics. 

In Figure |2| a comparison between the full dispersion relation m and the 5-th order Taylor 
expansion (jlip is shown for the physically realistic nonlocality = 200 [17|. It can be seen that (HU 
is a good approximation to the full dispersion relation in the limit of low wavenumber, as expected. 
However, due to the large factor in front of the term in the approximate dispersion relation HU 
the low wavenumber expansion rapidly deviates from the exact dispersion relation as k increases. 
Nevertheless, it qualitatively captures the key feature of the full dispersion relation, its non-convexity, 
so can be used for qualitative predictions of the effects of nonlocality on the nematic DSW behaviour. 
It is further seen from Figure [2] that the full phase velocity co/k is not monotone and has a minimum, 
which is also qualitatively captured by the long wave dispersion relation HU- The corresponding 
nonlinear equation with this linear dispersion relation, the fifth order KdV equation, will be derived 
in the next section. 

Let us now look at the opposite, dispersionless limit of the nematic system (I6|)-(l8|), which is 
described by the hyperbolic system of shallow water type [T] 


dp d , . 

= 0, 

(13) 

dv dv d9 


(14) 

dz dx dx 

= 0, 

9 

P 

(15) 
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These equations can be set in the Riemann invariant form 


V + 

V — 


2V2 

2V2 

7^ 


Vp 


V~p 


= constant 

= constant 


on 


on 


^ dx ,, _ 

dx ^ 


( 16 ) 

(17) 


The rarefaction wave seen in Fig. [T]can then be described by a centred simple wave solution of equations 
(1161) and (1171) in which the right-going Riemann invariant is constant. This solution will be presented 
in Section 05.11 


4 Fifth order KdV equation 

It has been shown that the nematic system ([1]) and ([2]) reduces to the Korteweg-de Vries (KdV) 
equation in the limit of small deviations from a background level [53( I22j . However, the physically 
large value of the nonlocality v m and the linked resonant wavetrain have major effects on the 
asymptotic analysis, which were not considered in this previous work. 

Indeed, assuming u ^ 1, but = 0(1), one has to retain the fifth order terms in the dispersion 
relation expansion (lllj) . which implies the necessity of keeping the fifth order dispersion term in the 
corresponding asymptotic KdV equation. The asymptotic reduction of the nematic equations to the 
KdV equation in the limit of small deviations from a background level uq will then be revisited, taking 
account of the large value of the nonlocality v. This asymptotic KdV equation will be derived from 
the hydrodynamic form of the nematic equations ([6])-([8|). Let us expand the hydrodynamic variables 
as 


Po + e^Pi {(, v) -h e'^P2 {C,ri) + ..., 

(18) 

rj) + e‘^V2(?, 0) + 0) + • • •, 

(19) 

— + e^ 6 'i(C, p) + p) + e® 03 (C, ??) + ..., 

(20) 


where po = 0 < e <C 1 is a measure of the deviation from the background uq, here = U 2 — U 1 . 

Also, ^ = e{x — Uz) and p = e^z are the usual stretched variables used to derive the KdV equation 
[I]. We also assume that all corrections to the equilibrium state p = po, v = 0, 9 = po/q vanish as 
\ C \ 00. 

Substituting the expansions ([TH]) and (pH]) into the director equation ([2]), we obtain at O(e^) 


9i 


Q 


( 21 ) 


and at 0(e'^) 


d^ 9 i P2 d‘^62 
2q 9^2 q 2q 8^“^ 


( 22 ) 


The term ve^ 62 ^^j 2 q is formally O(e^) and should appear at next order in the expression for @ 3 , as in 
[53]. However, this implicitly assumes that v = 0(1), which is not the case for experimental values of 
v. Hence, this term will be retained at O(e^). Treating ve^92i^(^l2q as a correction, equation (I22p can 
be solved for 02 to give 


2q 8^ q 8^^ 2q^ 8^"^ 


(23) 


Note that the last term in (j23p has to be retained as (j22[) implies that P 2 can be of O(z^), making the 
last term 0{u^e^). 
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Substituting the expansions (fT8l) - (l20]) into the “mass” and “momentum” equations Q and ([7]), 
we have at O(e^) 


dVi U dPi , ^dVi 2d Pi 


(24) 


respectively, on using ([2T]) for 6i. Compatibility between these two equations for Vi and Pi then gives 
the coordinate velocity U as 



(25) 


q 


Identifying Uq = p from Section [3l we see that U = c from the long wave expansion (llljl of the linear 
dispersion relation. 

Similarly, at O(e^) the mass and momentum equations ([6]) and ([7]) give 


dV2 


+ 


dPi 

dr] 


+ Ii 


dPi 


+ Pi^=0 

dC 


(26) 


and 


d^ d^ drj ^ d^ 


1 d^Pi 
4po d^^ 


= 0 , 


(27) 


respectively. 

It was shown in [221 [53] that substituting the leading order part of (1231) (the terms in brackets) 
into (PT)) and combining it with and ([75]) leads to the KdV equation. We now need to extend 
this derivation by including the higher order terms of (1231) . The problem we encounter is with the 
computation of the last term in (I23p as the correction P 2 cannot be computed separately at order O(e^), 
leading to equations (1261) and ([27)) . and a higher order approximation is required. This difficulty can 
be circumvented by suggesting a suitable ansatz for P 2 which is compatible with (1261) and ((271) . Let 


P2 = air 


d‘^61 

de 


vd'^Pi 
q 0 ^^ 


where a is a constant. Then substituting (l23l) and ([28]) into ([27]) we obtain, on using (1241) . 


(28) 


dV2 _ 1 \dPi 2U dPi auUd^Pi 

di po dp ^ po ^ di q d^^ 


(29) 


Substituting (1751) . ([75]) and (1771) into (1771) we obtain the fifth order KdV equation for Pi 


dPi 3 dPi U / ir q \ d^Pi iz'^e^pp d^Pi _ 

dp qU ^ d^ 4 Vg 4po^ d^^ Aq^U ^ d^^ 


(30) 


For the 5th order KdV equation (1301) to be consistent with the long wave expansion (|lip of the linear 
dispersion relation [T] we have to choose a = —1/4 (note that due to the scaling for ^ and p one has 
to replace {io — kc) —)• e^cu, k ^ ek in (1TT|) to make the comparison). We note that if the substitution 
([78l) were not compatible with equations (1761) and ([77]) . it would not be possible to obtain agreement 
for both dispersive terms in (1301) with the expansion of the nematic dispersion relation (llip using the 
single fitting parameter a. 

The 5th order KdV equation ([371) differs from that found in (53] [77] due to the Pi^^^^^ term, which 
arises at this order as v is large. The polarity of the solitary wave solution of the 5th order KdV 
equation (I3np depends on the sign of the coefficient of the Pi^^^ term. It is then clear that in the 
nonlocal regime with u large the solitary wave solution of the defocusing nematic equations ([T]) and 
([ 2 ]) is a bright solitary wave, rising above a background level, rather than the usual dark solitary wave 
of the defocusing NTS equation, which the nematic equations become in the limit z/ —)• 0. 
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Figure 3: Numerical solution of 5tlr order KdV equation (l3T]) for wq = 0.5 and 7 = 0.05. red (dashed) 
line: initial condition at t = 0 , blue (solid) line: solution at t = 20 . 


Although the fifth order KdV equation (I30|) has a limited range of validity as an asymptotic, 
quantitative model for nematic DSWs, it provides major qualitative insight into their dynamics by 
capturing the effect of resonant radiation. To illustrate this, we solved numerically the normalised 5th 
order KdV equation 


dw dw d^w d^w 


(31) 


for sufficiently small 7 > 0. Equation (|31ll has been derived in several physical contexts, including 
magnetoacoustic waves and capillary-gravity waves of small amplitude when the Bond number is close 
to, but just less than, 1/3 (see e.g. [l2] and references therein). Radiating solitary waves solutions of 
()3ip were discovered by Kawahara [54] and then studied analytically and numerically in a number of 
papers (see e.g. [ISl lUllSO] and references therein). 

Let us consider the 5th order KdV equation (j31|) with the initial condition to = 0, x > 0 and 
w = Wq, X < 0, so that a DSW is generated. Due to the non-convexity of the dispersion relation for 
dsn), there is the possibility of energy exchange between long and short waves propagating with the 
same phase velocity (see Figure [2|), so this DSW is expected to generate a resonant linear wavetrain 
propagating ahead of it |5T]. Such a radiating KdV DSW is displayed in Figure [3l The solution shown 
in this figure has strong similarities to the radiating nematic DSW solution of Figure [T](a). However, 
the resonant wavetrain of the nematic solution is more uniform, which is due to the smoothing effect 
of the large nonlocality i'. 

In conclusion, we note that, although the known theory of radiating solitary waves provides some 
intuition as to the counterpart radiating DSW solution, the major contrasting feature of radiating 
DSWs is the fact that the lead solitary wave of the radiating DSW remains steady, while an isolated 
radiating solitary wave is intrinsically unsteady due to the radiation carrying away the solitary wave’s 
energy [48] . 
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5 Dam break problem for the nematic system 

The solution of the Riemann problem ([I])-(l4]) for the nematic system generically consists of three 
distinct parts: a rarefaction wave, a (bright) DSW and a radiative wavetrain (see Figure (Ha)). The 
rarefaction wave was analysed in m, so below we only briefly outline the relevant results. Our main 
attention in this section will be on the DSW on the intermediate level U 2 and the resonant wavetrain 
generated by it. 


5.1 Rarefaction wave 


The solution displayed in Figure (Ha) shows that there is an expansion wave linking the initial level 
tt 3 behind the DSW to the level U 2 on which the KdV DSW sits. This expansion wave solution has 
already been determined [22], so only the relevant details will be given here. 

The expansion wave linking the initial level behind to the intermediate level U 2 = can be 
found as a simple wave solution of the Riemann invariant equations (1161) and (HID as [22] 

X - V 2 u^ 

2 ^ Vq 

(32) 

V 2 


V~P = 


3 V 2 . 


U3, 

2V2n,3 _ ^ 


^(2«3-3yft) <f <s+ 


and 


V = 


0 , 

2-\/2tt,3 I 2x 

^ («3 - s/m) 


- < — 
z 


\/ 2 us 


^/9 


-^ < f < ^ {2 u3 - 3Vp5) 

^(2«s-3VS)<i<s+ 


(33) 


Here s+ is the velocity of the lead soliton of the KdV DSW, which lies at the leading edge of the 
intermediate shelf. The simple wave solution (|32p linking the initial level U 3 and the intermediate 
shelf U 2 will be used for the comparisons with numerical solutions in Section (6] 

The velocity of the leading edge of the DSW on the shelf U 2 needs to be determined. In contrast 
to the KdV and defocusing NLS equations, this velocity is not determined by the conservation of the 
Riemann invariant on C_ across the DSW [28], but by the classical shock jump condition. We note 
here that the occurrence of the classical shock conditions in a conservative dispersive hydrodynamics 
was observed earlier in numerical simulations of large amplitude shallow water DSWs m and optical 
DSWs in photorefractive media [56] (see also m) and, very recently, in the context of radiating 
dispersive shock waves governed by the defocusing NLS equation modified by third order dispersion 
[HKniiadT]. This remarkable generic phenomenon requires further analytical study. 

The non-dispersive equations (fT3l) - (fT^ have the jump conditions [I] 

i„2 I 2p2 

P 2 V 2 , 2^2 + q 

s+ =- and = — 


2pi 

g 


P2 — Pi V2 

as ahead of the shock, p = pi and u = 0 and behind the shock, p = P 2 and v = V 2 - 
between these equations gives 

2 p2- Pi , 2 p2 


V 2 = 


and = 


(34) 

Eliminating 

(35) 


The expansion fan solution (1331) also gives an expression for V 2 in terms of the intermediate level U 2 - 
Matching this and (1351) then gives that this intermediate level U 2 = P 2 is the solution of 

U 2 — + 2 + 2u‘f^ u\ — ^u^u\u 2 — u\ + 2u^u\ = 0 (36) 

with ui < U 2 < U 3 . For the particular case ui = 0, ri 2 = (2 — \/2)u3. Also, as ui u^, U 2 —>■ 
{u 3 + ui) /2, which is the value obtained by conservation of the Riemann invariant (I17p on C- 
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5.2 Dispersive shock wave: lead solitary wave 

In previous work [22] the DSW solution of the KdV equation |25l [26| was used for the DSW on the 
intermediate shelf U 2 - While this was found to give good agreement with numerical solutions for 
values of ui near U 3 , significant disagreement was found for values of ui away from U 3 . As discussed 
above, this is due to the velocity s+ of the front of the full nematic DSW not being well determined 
by the velocity of the leading edge of the asymptotic KdV DSW. The reason for this behaviour is 
that the DSW is subject to radiative losses due to the resonance with the co-propagating linear short 
wavelength waves, resulting in a rapidly oscillating wavetrain shed ahead of the DSW. For small initial 
steps the radiating DSW is described in the framework of the 5th order KdV equation (I30p . However, 
for general jumps the full nematic system should be used due to the 5th order KdV equation not being 
accurate in capturing large wavenumber dispersive behaviour (see Fig. [2|). 

We now need to relate the shock velocity s+ to the amplitude of the lead soliton of the DSW. 
Since the solitary wave solution of the full nematic system is not available, as an approximation we 
shall use the soliton solution of the standard KdV equation, that is (fSO]) without the 5th derivative 
term. On noting the scalings in the expansions (fT8]l - ([20]) and equating given by (1351) to the lead 
soliton velocity, this gives 

as = e A= -m- (37) 

V ^2 + 

The lead soliton of the KdV DSW itself is given by [22] 

|tt| = y/p = uo + e^Asech^/3(x — s+z) + ..., (38) 

where 

^- 

V2{2qyyay/V 

These results will be used in the next section to find a solution for the resonant wavetrain leading the 
KdV DSW seen in Figure (Ha) on identifying uq = ui. 


U 


qU 


Aq I6u1v 




(39) 


5.3 Resonant wavetrain 


Let us now consider the wavetrain ahead of the DSW. This wavetrain is generated due to a resonance 
between the long wave oscillations in the DSW and co-propagating short wavelength waves, as implied 
by the non-convexity of the linear dispersion relation (HOI) and discussed in Section 01 

In determining the structure of the resonant wavetrain we refer to Figure UKa) in which one can 
observe three regions of distinctly different behaviour: region (i) provides a transition from the lead 
soliton of the DSW to region (ii) which contains the (almost) uniform extended middle part of the 
wavetrain; and the front region (Hi) which brings the wavetrain down to the constant level u = ui 
and 6 = 0i, where 0i = |uip/g, see dlj). 

We start with the middle region (ii) for which the director 6 is close to 0i, 0 — 0i = and 

the wavenumber k is 0(1). Hence, the asymptotic dispersion relation (II2p applies and the dispersion 
relation for the resonant wavetrain is 

= 1/^2 + 201 ( 40 ) 


as the resonant wavetrain is on the background carrier wave uiexp(—2iQiz). Furthermore, the res¬ 
onant wavetrain in the region (ii) is then asymptotically described by the linear equation following 
from (HI) on setting 0 = 0i, 


,du 1 d^u 
* 2 5x2 


201U = 0. 


( 41 ) 
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We assume that the main resonance is with the lead soliton of the DSW, which we approximate 
by the KdV soliton (f38l) . Matching the phase velocity to the lead KdV soliton velocity (|35]) . we have 


(42) 


Cr = ^ 

which can be solved to give the wavenumber of the resonant wavetrain as 


k = kr = s+ + 


*+ - ^“1 


4 


The front of the resonant wavetrain moves at the group velocity Cg [T], which is 

Cg = (jj'r{kr) = K- 


(43) 


(44) 


These expressions for the asymptotic wavenumber of the resonant wavetrain away from the DSW and 
the velocity of the front of the wavetrain will be used in the solution for this wavetrain. 

The wavenumber (l43]l is real if ui < uic, where uic is the solution of 

2 Un 2 

'®+ ~ r, ~ —pUic 


(45) 


For ui above uic there is only a transient wavetrain ahead of the DSW [I]. This existence of a critical 
ui above which there is no resonant wavetrain is in agreement with previous work [22] in which the 
critical value was found to be uic = For tts = 1, g = 2 and v = 200 numerical solutions 

give the critical value uic = 0.69 [22]. For these parameter values the new modulation value (1451) 
uic = 0.648 is slightly below the numerical cut-off, while the previous modulation value uic = l/\/2 
is slightly above. It should be noted that numerical solutions do not show a sharp transition to no 
resonant wavetrain as given by (I45p . but a rapid transition from an upstream uniform wavetrain to 
none over a ui range of about 0.1. 

Above the critical value (|45h the resonant wavetrain ceases to exist. The DSW on the intermediate 
level U 2 then becomes the standard KdV type DSW and the approximate solution of [22] holds. The 
amplitude and velocity of the lead soliton of the DSW are, then cf. (l3^ . (f37l) . 


as = e A = uz-ui, s+ = \l-U 3 . 


(46) 


As equation ()4ip is linear, it does not allow the determination of the resonant wavetrain ampli¬ 
tude. For that, one needs to go beyond the approximation 0 = 0i in the wavetrain and include the 
(signihcant) variations of the director in the transition region (i) between the DSW and the uniform 
wavetrain region. 

As the phase velocity of the resonant wavetrain is the same as the (classical shock) velocity s+ (1351) 
of the lead soliton of the DSW, to determine the solution for the wavetrain in the transition region 
we will use the moving coordinate C, = x — s+z. By inspection of the numerical solution of Figure 
dDa) it is reasonable to assume that the approximate director solution in the transition region is given 
by the lead solitary wave of the DSW, so that from equations and ([SSI) of the KdV expansion of 
Sections 0] and 55.21 we have 

0 = — -I- e^^^Asech^/3C- (47) 

q q 

The ansatz (I47p transforms the equation for the electric field ([T]) into a linear, variable coefficient 
equation whose solution can be sought in the form 




( 48 ) 
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where (t(C) and Ur{C,z) are the phase correction due to the variable coefficient and the wavetrain 
amplitude, respectively. To be consistent with the director (I47h Ur = O(e^) as it is proportional to 
the jump height U 2 — ui. Substituting (H71) and (H8|) into the electric field equation ([T]) we have 


dUr dUr 1 d‘^Ur ( 4e^Ui 

+ ct'uis+ - Asech^/?C — 0- 

2 q 


Asech^/3C — s+a' — -a” + -a'~ ) Ur 


(49) 


We now choose the phase correction (t{C) so that the relation 


s+a' — -a'‘^ = —— Asech.'^/3C (50) 

z q 

is satisfied. Then, on using (I50p to leading order in e, we obtain from (1491) the equation for the variation 
of the wavetrain amplitude in the transition region as 


i 


dur 

dz 


i{s+ - cr')-^ + 


dUr . 1 8“^ Ur 


dC 2 dC^ 


= 0 . 


(51) 


In deriving this equation, we have noted that a" is higher order in e ( since (5 ~ see (l39]l ). 

Using the numerical solution (see Fig. UKa)) and the soliton solution (p8]l as a guide to the structure 
of the transition region we shall look for the solution of equation (1511) for Ur as fast (scaled as 0(1)) 
oscillations with a slowly varying (scaled as 0(/3“^)) envelope. We then seek a WKB solution of the 
form 

Ur = W{X, ^ (52) 

where the slow variables are X = PC and Z = Pz. This WKB expansion is valid if l/\/T <C e ^ ^/u, 
which holds as u is large. The first inequality is due to using the hrst two terms of the KdV expansion 
(f20l) for the director (ITFl) and the second is required for the validity of the WKB form (l5^ . Substituting 
the WKB form (j52|) into equation (j51|) gives the eikonal equation 


and the transport equation 


dpj 


dijj 


1 [ dip Y / l\--r 


dW (dip . 


dW 


1 d‘^ip 

2 ^ 


W. 


(53) 

(54) 


We note that the group and phase velocity argument gave that as the resonant wavetrain ap¬ 
proaches the wavefront at x = CgZ, it becomes a uniform wavetrain of wavenumber kr and frequency 
kr/2 + 201 [22]. We then find that the solution of the eikonal equation ([531) is 


/I \ A) ^\/ 

iP = krX- [-kl - s^k) Z - —^^-^^tanhX, ~P = ^.P= 

We note that the phase correction (f5^ becomes infinite as kr ^ Sj^. This is expected as the group 
velocity of the front of the resonant wavetrain is kr- When the velocity of the lead soliton of the KdV 
DSW is greater than the group velocity, the wavetrain cannot propagate away from the DSW. There 
is then no upstream resonant wavetrain, with only a small amplitude transient being present [22| . 

To solve the transport equation (|54[) the resonant wavetrain leading the KdV DSW must be 
matched to the intermediate shelf, so that W = Wq = = U 2 — ui at V = 0 on noting the full 
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solution (I48p for u. Then using the eikonal equation solution (I55p . the solution of the transport 
equation (jMP is 


W = Wo 


1 + 


2uie^krA 


rsech^X 


1 + 


2uie^kj.A 


1 -1 


qs+{kr - s+)2 J [ qs+{kr - s+)2 

The height of the resonant wavetrain exponentially approaches the constant value 

. 2 ; . 


Wc = Wo 


1 + 


2uie^krA 
qs+{kr - s +)2 


(56) 


(57) 


as the front of the wavetrain at x = CgZ is approached, so that the total height of the envelope of the 
resonant wavetrain in the region (ii) is given by 


CL'p 


ui + Wo 


1 + 


2uie^krA 
qs+{kr - s +)2 


1 -1 


(58) 


Finally, we describe region (iii) of the resonant wavetrain which brings it down to the initial level 
ui (see Figure dK a)). In the region of this front, as for the uniform middle region (ii), we approximate 9 
by 0 = 01 = |uip/g, so that the linear equation (jUp holds. If we use a moving coordinate Qg = x — CgZ 
moving with the velocity of the front, the electric field is governed by 


du du 1 d‘^u 

'Tz - + 2 ^ - = "■ 

To match with the initial level ahead, we seek a solution of the form 


(59) 


(60) 


so that Uf IS the solution of 


.duf . duf ^ I d'^Uf ^ 


dz 


dCg ' 2 dQ 


(61) 


To match with the uniform wavetrain behind, we have the boundary condition \uf \ = a^ — ui at Cg = 0. 
The linear equation (| 6 ip can be solved using Laplace transforms to give the Fresnel integral solution 


Uf 


2{ar Ui) 


POO 

y/2z 


dt. 


(62) 


6 Comparison with numerical solutions 

In this section, full numerical solutions of the nematic equations ([T|) and ([ 2 |) will be compared with the 
modulation theory solutions of Sections 1^ 15.1115.21 and 15.31 The numerical solution of the electric held 
equation ([T|), which is of NLS-type, was obtained using the pseudo-spectral method of Fornberg and 
Whitham [26], modihed to improve its accuracy and stability [58|, but without the boundary damper 
due to the non-zero boundary conditions. These improvements include using a 4th order Runge-Kutta 
scheme to propagate forward in z, resulting in higher accuracy, in Fourier space, rather than in real 
space, resulting in improved stability |SS|. The numerical solution of the linear director equation 
([^ was obtained using a spectral method [59]. This numerical scheme is discussed in [60]. For the 
numerical solutions of this work 32768 points were used for the FFT with a x domain of length 8192.0 
and a z step of Az = 0 . 002 . The x domain was chosen long enough so that the waves at the numerical 
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Figure 4: Comparison between numerical solution of nematic equations ([I]) and ([ 2 ]) and the modulation 
theory solution of Sections 135.11135.21 and 135.31 for us = 1.0, ui = 0.5, q = 2 and = 200. Initial 
condition for |rt| {z = 0): red (short dashed) line; initial condition for 9 [z = 0): green (long dash) 
line; numerical solution for |rt| aX z = 300: blue (solid) line; numerical solution for 9 at z = 300: pink 
(dotted) line; modulation solution: black (dot-dot-dash) line. Only the lead soliton of the modulation 
theory DSW solution is shown. 


boundaries generated by periodicity were far from the region of interest. Finally, the initial condition 
([3]) was smoothed using tanhx/kF to avoid spurious numerical effects due to large x derivatives, with 
W = 1 found to be suitable. 

Figure H] shows a comparison between the numerical solution of the nematic equations ([T]) and 
(|^ for ri 3 = 1.0 and ui = 0.5 at z = 300 for 5 = 2 and v = 200. For clarity, in these figures 
only the upper envelope of the resonant wavetrain ()48]) . (1521) and ([56]l and the upper envelope of the 
Fresnel front (I62p are shown. It can be seen that there is very good agreement in general between the 
numerical solution for the electric field |u| and the modulation theory solution of Sections 35.11 35.21 
and 35.31 In particular, there is excellent agreement for the position of the lead soliton of the DSW, 
which is the same as that of the trailing edge of the resonant wavetrain. This is in contrast to the 
result of previous work | 22 ] in which this position was determined by the velocity s+ = y^ 2 /gu 3 of the 
lead soliton of the standard KdV DSW solution [25l[26], resulting in the DSW leading edge being at 
X = 300 for the parameters of Figure [H noting that the numerical position is x = 247.5. It is then 
clear that the shock velocity (I35p determined from the shock jump conditions for the non-dispersive 
“shallow water” equations (fT^ - (fT5|) and giving x = 265.9 for the lead soliton at z = 300 yields much 
better agreement with the numerical solution for the position of the leading edge of the DSW than 
the velocity determined by the KdV DSW solution. The differing length scales of the KdV DSW 
{0{^/l')) and the resonant wavetrain (0(1)) can be clearly seen. The major disagreement is that the 
front of the numerical resonant wavetrain has more structure than the linear Fresnel integral solution 
of Section 35.31 However, the Fresnel integral solution gives the correct spatial extent of the transient 
front of the resonant wavetrain. Furthermore, if the Fresnel integral solution is shifted so that it starts 
ahead of the rise in the numerical front, it is in very good agreement with the numerical front. 

The other noticeable disagreement between the numerical and analytical solutions is the amplitude 
of the lead soliton of the DSW. The amplitude of the DSW in the electric held u is generally under¬ 
predicted by the KdV theory of Section 35.21 which was based on the classical shock speed (|35]l . 
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Figure 5: (a) Height of lead soliton of the DSW. Numerical solution of ([T]) and m- red pluses; 
analytical solution Ug = ui + e^A (f37)l and (ll 6 l) : green (dashed) line, (b) Height U 2 of intermediate 
shelf. Numerical solution of ([T]) and m- red pluses; modulation solution (l36]l : green (dashed) line, 
(c) Height ttr of resonant wavetrain as a function of the upstream level ui. Numerical solution of ([T|) 
and ([2]): red pluses; WKB solution ([58]): green (dashed) line, (d) Comparison for leading and trailing 
edges of resonant wavetrain at z = 300. Numerical trailing edge: red pluses; trailing edge x_ = s+z 
given by the classical shock speed (H 2 ]) : red (solid) line; numerical leading edge: green crosses; 
leading edge = CgZ defined by the group velocity (jUj): green (dashed) line; trailing edge of [ 22 ] 
given by the soliton speed ([l 6 ]l in the standard modulation solution for the KdV DSW: blue (dotted) 
line. The other parameter values are 7 x 3 = 1, q = 2 and 1 / = 200. 


However, this approximation yields good agreement for the DSW in the director 6, given in the KdV 
approximation by Eq. ()47l) . This is in contrast to the results of |22] for which the standard DSW 
solution of the KdV equation was used to determine the DSW on the intermediate level U 2 - The 
results of m strongly over-predicted the height of the bore in the director, this major discrepancy 
being fixed in the present theory. Finally, it can be seen that under the resonant wavetrain there is a 
slight rise in the director above 0 = 0i due to 0(zx“^) corrections in the asymptotic expansions. These 
higher order corrections will be dealt with in future work based on a full description of a resonantly 
radiating DSW. 

The agreement between the modulation theory and numerical solutions is further quantified in 
Figure [5j Figure [5Ka) shows a comparison of the height (background plus amplitude) of the lead 
soliton of the DSW as given by numerical solutions and the modulation solution (a^ = e^A + ui), 
using (1371) for the amplitude below the cut-off (1451) and (HBI) above. The choice of the total height 
rather than amplitude for the comparisons is due to the soliton background being not clearly defined 
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in the numerical solutions (see Fig. UKa)). Furthermore, the amplitude of the lead wave oscillates 
slightly due to its interaction with the resonant wavetrain, so the figure shows the average amplitude. 
The numerical solution clearly shows the predicted different DSW behaviours above and below the 
resonant wave cut-off, which was not predicted in [ 22 ] for which the height was the constant value 
(1461) for the whole range of ui. Thus, one can see that the KdV soliton height based on the classical 
shock wave speed is in broad agreement with the numerical values. The appropriateness of using the 
classical shock wave velocity to determine the intermediate shelf height (I36p is quantified in Figure 
[5Kb). It can be seen that (I36p is in excellent agreement with the numerical height, except for a slight 
discrepancy as ui —)• 0. This is due to the intermediate shelf disappearing as the dam break solution 
for ui = 0 is approached [ 22 ] . 

Figure [5](c) shows a comparison between the height of the resonant wavetrain obtained numerically 
and the modulation solution height ()58p . There is excellent agreement between these heights, except 
towards the cut-off near ui = 0.7. This is due to the discrepancy between the numerically found 
cut-off and the modulation theory prediction. 

Finally, Figure[5Kd) shows a comparison for the leading and trailing edges of the resonant wavetrain. 
It can be seen that there is excellent agreement for the position of the trailing edge, even up to the 
cut-off. Previous work m predicted the constant (i.e. independent of ui) velocity (1461) for the trailing 
edge which was defined by the value U 3 alone. Figure [5Kd) clearly shows that the present theory 
based on the classical shock speed for the leading edge of the DSW is in much better agreement. The 
agreement for the leading edge is reasonable above ui = 0.5 as the cut-off is approached, but is poor 
as ui decreases. There are a number of reasons for this. The position of the trailing edge is clearly 
defined by the peak of the lead soliton of the KdV DSW. While the theory of Section 55.31 predicts a 
precise location for the leading edge of the resonant wavetrain, it can be seen from Figs, dja) and d] 
that this is not the case for the numerical solution. There is no clean boundary between the resonant 
wavetrain and its front. There is an extended transition between the two. For the comparison of Figure 
[SKd) the start of the hump before the monotonic decrease of the front was chosen as the leading edge 
position of the wavetrain. Finally, the present modulation solution underpredicts the cut-off point 
for the resonant wavetrain (at ui = 0.648 compared with the numerical value 0.69 for the parameter 
values of Fig. [5]), leading to the disagreement as the cut-off is approached. 

Part of the reasons for these discrepancies is due to the major simplifying assumption adopted 
in the present theory, according to which the generation of the wavetrain is dominated by a single 
resonance with the lead soliton of the DSW. A more advanced theory including internal resonances 
with the other components of the DSW is needed to achieve better agreement with numerical solutions. 

There is one more feature of the resonant wavetrain which complicates its analysis for large initial 
jumps. As the initial level ahead ui decreases the electric field u eventually vanishes at a point, 
termed the vacuum point m- For sufficiently large initial jumps the vacuum point occurs within the 
resonant wavetrain, so that the lower envelope becomes non-monotone. It was shown in |37j that for 
the defocusing NLS DSW there is a singularity in the phase v at the vacuum point itself. Although 
the resonant wavetrain for the nematic system ([T|) and (( 2 |) is asymptotically described by the linear 
equation (1411) rather than the defocusing NLS equation, numerical simulations show that the vacuum 
point in the wavetrain has qualitatively similar properties to the vacuum point arising in the large 
amplitude NLS DSW [37] . In particular, such a DSW has a non-monotone lower envelope (see Figure 
[ 6 Ka)) and exhibits a phase singularity at the vacuum point, see Figure [ 6 Kb). The WKB solution of 
Section 55.31 gives that the lower envelope of the resonant wavetrain has height 


ai = ui-Wc = ui- Wo 


1 + 


2uie^krA 
qs+{kr - s +)2 


1 -1 


(63) 


For U 3 = 1.0, u = 200 and q = 2, it is found that ai vanishes when ui = 0.2416. Numerical solutions 
of the nematic equations ([T|) and ([ 2 |) show that for these parameter values a vacuum point first occurs 
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Figure 6 : Numerical solution of the nematic equations m and (I 2 |) for = 1 . 0 , Ui = 0 . 1 , q = 2 and 
u = 200. (a) Initial condition for |ri| (z = 0): red (short dashed) line; initial condition for 6 (z = 0): 
green (long dash) line; numerical solution for |u| at 2 : = 300: blue (solid) line; numerical solution for 

6 at z = 300: pink (dotted) line, (b) phase v at z = 300. 

when ui = 0.22. A full analysis of the solution after the vacuum point is reached is beyond the scope 
of this paper. Full Whitham modulation equations would be required for a proper analysis after the 
vacuum point m- 

7 Conclusions 

The Riemann problem for the equations governing the propagation of a coherent optical beam in a 
defocusing nematic liquid crystal has been studied. It was found that in the highly nonlocal limit 
the DSW, which comprises a major part of the Riemann problem solution, is drastically different 
to the DSW solution of the defocusing NTS equation, to which the nematic equations reduce in the 
small nonlocality limit, that is —)• 0. There are two major differences: (i) the nematic DSW is of 
positive polarity with a bright soliton at its leading edge; (ii) it is preceded by a short wavelength 
resonant wavetrain. To clarify this structure, it was shown that in the limit of small deviations from 
a background, the nematic equations reduce to a KdV equation with a fifth order derivative, the 
Kawahara equation. This fifth order KdV equation is known to have a resonance between its solitary 
wave solution and linear radiation. The present work shows that this resonance extends to a resonance 
between the DSW and linear radiation. A modulation theory was developed to derive solutions for 
the resonant wavetrain and its front. In contrast to previous work [22], it was found that the leading 
edge of the DSW was determined by the classical shock jump condition, which is non-standard for 
DSWs [28|. Excellent agreement was found between the major part of the modulation theory solution 
and full numerical solutions of the nematic equations. However, there are some discrepancies. Part of 
the observed discrepancies can be addressed by applying a more complete modulation theory for the 
DSW description. When the higher order dispersion term is a small perturbation, as in the Kawahara 
equation (f3T]l with 7 <C I, Whitham theory for perturbed integrable equations [HI [62| provides an 
appropriate analytical framework for the description of the DSW evolution. 

The present work leaves open a number of issues. As already discussed, resonances between DSWs 
and radiation is an issue which has received little attention to date with the theory and solution 
methods only starting to be developed [snini, in contrast to the corresponding resonant interaction 
between solitary waves and radiation (see e.g. [ISl HHl [50l [SU 03]). As the nematic equations are 
generic and similar equations arise in a number of fields, this resonant interaction deserves in depth 
treatment. A proper analysis of possible resonances between DSWs and radiation is an open question 
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which deserves further treatment. This wiii be the subject of further work. 
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